Nonlinear interaction effects in a strongly driven optomechanical cavity 
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We consider how nonlinear interaction effects can manifest themselves and even be enhanced in 
a strongly driven optomechanical system. Using a Keldysh Green's function approach, we calculate 
modifications to the cavity density of states due to both linear and nonlinear optomechanical in- 
teractions, showing that strong modifications can arise even for a weak nonlinear interaction. We 
show how this quantity can be directly probed in an optomechanically-induced transparency type 
experiment. We also show how the enhanced interaction can lead to nonclassical behaviour, as 
evidenced by the behaviour of (72 correlation functions. 

PACS numbers: 



Introduction- The field of cavity optomechanics in- 
volves understanding and exploiting the quantum in- 
teraction between a mechanical resonator and photons 
in a driven electromagnetic cavity. It holds immense 
promise for both fundamental studies of large-scale quan- 
tum phenomena as well as applications to quantum infor- 
mation processing and ultra-sensitive detection, and has 
seen remarkable progress in the past five years. High- 
lights include the use of radiation pressure forces to cool 
a mechanical resonator to close to its motional ground 
state [II [2] and experiments where the mechanical motion 
causes squeezing of the light leaving the cavity |3l 1^ . 

As remarkable as this progress has been, it has re- 
lied on strongly driving the optomechanical cavity to 
enhance the basic dispersive coupling between photons 
and mechanical position. While the resulting interac- 
tion can be made larger than even the dissipative rates 
in the system [SHZ], it is purely bilinear in photon and 
phonon operators. As a result, it cannot convert Gaus- 
sian state inputs into non-classical states or give rise 
to true photon-photon interactions. Recent theoretical 
work has addressed physics of the nonlinear interaction 
in weakly driven systems [SI [5] ■ Unfortunately, one finds 
that effects are suppressed by the small parameter g/uJM- 

In this paper, we now consider nonlinear interac- 
tion effects in an optomechanical system that (unlike 
Refs. [3 [9]) is also subject to a strong laser drive; we con- 
sider effects of this driving beyond simple linear-response. 
We find somewhat surprisingly that one can use the 
strong drive to enhance the underlying single-photon in- 
teraction. Using non-equilibrium many-body perturba- 
tion theory (based on the Keldysh technique (see, e.g., 
[4])), we calculate how these effects modify the cavity 
density of states, and hence the cavity's response to an 
additional weak probe laser. This response is exactly 
the quantity measured in so-called optomechanically- 
induced transparency (OMIT) experiments [HI fTTMl3j . 
We find striking modifications of the OMIT spectrum, 
effects which can be attributed to the nonlinear interac- 
tion causing a hybridization between one and two polari- 



ton states (with the polaritons being joint mechanical- 
photonic excitations). We also find the possibility of 
enhanced polariton-polariton interactions, which lead in 
turn to non-classical correlations (as measured by a §2 
correlation function). 

System- The standard Hamiltonian of a driven op- 
tomechanical cavity \s H — Hq + -ffdiss, with (h = 1) 
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(1) 
Here a is the cavity mode (frequency wc, damping rate 

k), 6 is the mechanical mode (frequency wm, damping 
rate 7), and g is the optomechanical coupling. iJdiss de- 
scribes dissipation of photons and phonons by indepen- 
dent baths; ain(i) is the amplitude of the drive laser. 

We consider the standard case of a continuous-wave 
drive (i.e. ain{t) oc e^*""-'), and work in a rotating frame 
at the laser frequency cjl- We further make a displace- 
ment transformation, writing a = e"**^""* la + dj, where 

a is the classical cavity amplitude induced by the laser 
drive. Letting A = lo-^ — loq, the coherent Hamiltonian 
now takes the form Hq = Hi + H2 with 

Hi = -Ad^d + ujMb^b + G{d + d^){b + b'<), (2) 
H2 - gSd{b + P). (3) 

G — ga is the drive-enhanced many-photon optomechan- 
ical coupling; we set g, a > without loss of generality. 

The most studied regime of optomechanics is where 
a ^ 1 and g ^ KjWm- It is then standard to neglect 
the effects of H2- In the absence of any driving, a sim- 
ple perturbative estimate suggests that the effects of H2 
enter as g^/wM, where the factor of ojm corresponds to a 
virtual state with one extra (or one less) phonon. This 
conclusion can be made more precise by exactly solving 
the coherent, undriven system using a polaron transfor- 
mation [HHH- Thus, in this standard regime, one can ig- 
nore H2 , leaving only Hi , which is easily diagonalized as 
Hi = Y^cr^^ Ecc\c„. Here c+,_ describe the two normal 
modes of the system. As these modes have both photon 



and phonon components, we refer to them as polaritons 
in what follows. Their energies are: 
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For A ~ —ujm and G > k, 7, the polariton energy split- 
ting can be resolved experimentally |5l-(7]. 

Polariton interactions- Unlike previous work, we wish 
to retain the effects of the nonlinear interaction H2 , but 
also consider the effects of a large drive (and hence a large 
many-photon coupling G). To proceed, we will treat the 
effects of H2 in perturbation theory. We use a Keldysh 
Green's function (GF) approach which is able to describe 
the non-equilibrium nature of the system. The linear 
Hamiltonian Hi along with the dissipativc terms in i^diss 
define the free GFs of the system, which describe the 
propagation of polaritons in the presence of dissipation. 
Written in the polariton basis, the nonlinear interaction 
H2 gives rise to number-non-conserving interactions. 
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where the coefficients g„lj,„n oc g [TT. Note normal 
ordering H2 in terms of polariton operators introduces 
small quadratic and linear terms which modify the diag- 
onafized form of Hi (see EPAPS for details [14]). 

We start by considering how single-particle properties 
are modified by the nonlinear interactions; such proper- 
ties can be directly probed by weakly driving the cavity 
with a second probe laser (i.e. an OMIT experiment fTTI - 
[13]) or by measuring the mechanical force susceptibility. 
Understanding these properties amounts to calculating 
the self-energy I][aj] of both the polaritons due to J?2- 
We have calculated all self-energy processes to second or- 
der in g. Our approach captures both the modification of 
spectral properties due to the interaction (i.e. the modifi- 
cation of the cavity and mechanical density of states), as 
well as modifications of the occupancies of the mechanics 
and cavity. While our approach is general, we will focus 
on the most interesting case of a high mechanical quality 
factor 7 ^ wmj a cavity in the resolved sideband regime 
wm > K, and a strong cavity drive, G > k. 

Our full second-order self-energy calculation finds that 
for most choices of parameters, the polariton self-energies 
scale as g^/uJu and thus have a negligible effect for 
the typical case where g ^ ujm- However, effects are 
much more pronounced if one adjusts parameters so that 
E^ = 2£'_. This condition makes the term in H2 which 
scatters a -|- polariton into two — polaritons (and vice- 
versa) resonant. It can be achieved for any laser detuning 
A in the range (— 2wm, — wm/2) by tuning the amplitude 
ciin of the driving laser so that the many-photon optome- 
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FIG. 1: Main: + polariton resonance in the cavity density of 
states, for various values of the nonlinear interaction strength 
g (as indicated), as obtained from Eq. ( 12 1 (with the inclusion 
of energy shifts from Hnr ^14J- For all plots, the laser drive is 
at the red sideband A = —ljm, and G = O.Sojm to ensure the 
resonance condition _E+ — 2E-; we also take lom/k, — 50, T = 
and 7 = 10~*«;. The peak splitting signals the hybridization 
of a -I- polariton with two — polaritons. The dashed curve is 
the result of a master-equation simulation for g = n [14j . 
Inset: Full density of states, same parameters, showing the 
asymmetry between -|- and — polariton resonances. 



chanical coupling G ~ Grcs, where 



G„s[A] = J4A4 - UA^ojIi + 4Ljti / lO^A^ . (6) 



In this regime, the dominant physics is well described 
by the approximation Hq ~ H^s with 
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The second term in H^s corresponds to making a 
rotating-wave approximation on the nonlinear interac- 
tion H2 in Eq. ([s]), retaining only the resonant process; 

g — g^ (_ ex (7 is the corresponding interaction strength 

(see inset of Fig.[2]to see how g varies with A). The terms 
in Hnr describe the small (i.e. ex g"^ juiyi) residual effects 
of the non-resonant interaction terms in Eq. ([Sl; we treat 
them via straightforward second-order perturbation the- 
ory (i.e. a Schrieffer- Wolff transformation) . They play no 
role in the extreme good-cavity limit ujyi ^ k [14j . 

Green functions for resonant nonlinear interactions- 
Focusing on the resonant interaction regime defined by 
Eq. (|6|, and using the simplified Hamiltonian in Eq. ([7|, 
we obtain simple expressions for the retarded GFs of the 
system. The retarded photon GF in the displaced, rotat- 
ing frame is defined as 



G 



ddV 



dt9{t)([d{t),d\0)])e 



(9) 



with similar definitions for tfie polariton retarded GF 
G^^H (a = ±). As usual, pdH = -Im G^^M/t^ de- 
scribes the cavity density of states; G^^[w] also deter- 
mines the reflection coefficient in an OMIT experiment 
(see Fig. [2]). A standard linear response calculation [T3] 
shows that the elastic OMIT reflection coefficient is given 
by rfajpr] = 1 — zKcpG^^[a;pr], where ujpi is the frequency 
of the weak probe beam, and Kcp is the contribution to 
the total cavity k from the coupling to the drive port. 

In the limit of interest where k <C Eg., there 
are no off-diagonal polariton GFs or self-energies [H] . 
As a result, G^^^[uj] will be given as G^^[w] = 

E, [c^Gtioj] + Dg [Gth^]] *) , where the change-of- 
basis coefficients Ca:Dg are given in [Hj. The Dyson 
equations for the polariton retarded GFs are 



G^.H 



Eg+inj2- 



s;vi^ 
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(10) 



where k„ is the effective damping rate of the a polariton 
|14j . Using the effective Hamiltonian in Eq. (l7|, a stan- 
dard Keldysh calculation yields that to second order in 
5, the polariton self-energies take the simple forms: 






2g2(l + 2n_) 

cj — 2E_ + iK^ ' 

4g^(n_ — n+) 
w - {E+ - E^) + i(K+ + K-)/2 



(11a) 
(lib) 



Here, n„ is the effective thermal occupancy of the a po- 

We have 



'rU-J\ ^ 



lariton [T3]; for g = 0, we have 
taken the limit g/uiM — >■ 0, and hence neglected the ef- 
fects of the non-resonant terms Hnr in Eqs. (11); the 



explicit corrections due to these terms are given in the 
supplemental information [14j . 



Eqs. ( 11 1 are central results of this work. Eq. ( 11a I de- 



scribes the fact that a single -f polariton can resonantly 
turn into two — polaritons, and describes the hybridiza- 
tion between these states that occurs for large enough 
g. To see this explicitly, we consider the case of exact 
resonance (i.e. E^ = 2E_) and write 



g5[ 



Ij: 



1 ■ 2k——k-\- 

1 - 'V iS^ + 



2^^^-E++^^-^^^^+r,5, 



(12) 



5+ = Y^2g2 (1 + 2n_) - (2k_ - k+f /16. (13) 

For g ^ K, we see that the + polariton GF has two 
poles, corresponding to the new hybridized eigenstates. 
We stress that these eigenstates do not correspond to a 
fixed excitation number. Note that unlike the undriven 
system [HI [3] , the effects of the nonlinear interaction can 
be significant even if g ^ wm- Also note that the reso- 
nant coupling between |-|-) and | ) states is enhanced 

at finite temperature by a standard stimulated emission 
factor (1 -|- 2n_). The form of this GF and self-energy 
are reminiscent to the photon GF for ordinary OMIT, 



where a photon can resonantly turn into a phonon [llj : 
however, that effect does not involve any temperature- 
dependent enhancement. S in Eq. (lib) describes a 

process where the propagating — polariton of interest in- 
teracts with an already-present — polariton to turn into 
a -f. As this process requires an existing density of po- 
laritons, it is strongly suppressed at low temperatures. 

We note that it is possible to use resonance to en- 
hance the nonlinear optomechanical interaction without 
strong driving, if one instead considers a system where 
two cavity modes interact with a single mechanical res- 
onator [TSHTT] . Our approach has the benefit of only re- 
quiring a single cavity mode; further, for drive detunings 
near A = — wm, it also has a natural resistance against 
mechanical heating, as mechanical contribution to the 
polariton temperature scales as "fnth/n, where fith is the 
mechanical thermal occupancy. While a low temperature 
is not essential for the density-of-states effects described 
above, it is essential for the correlation effects discussed 
below. Finally, for superconducting microwave cavities, 
the cavity linewidth k has a strong contribution from 
two-level fluctuators, and thus improves if one strongly 
drives the cavity (as in our scheme). 

Red-sideband drive- For a detuning A = — wm, the 
polariton resonance occurs when G — O.Swm- For this 
detuning, both polaritons are almost equal mixtures of 
photon and phonon operators. One finds Ko- = (k + 7)/2, 
and that the resonant interaction strength g ~ —0.37g. 
Because Hi does not conserve the number of photons and 
phonons, the polaritons are not eigenstates oid^d+b^b; as 
a result, even at zero temperature, the effective thermal 
occupancies scale as fi^ ex (G/wm)^ ^ 1 [H]- The inset 
of Fig. [l] shows the evolution of the cavity density of 
states for these parameters as g is increased from zero. 
For g = 0, one sees two symmetric peaks corresponding 
to the two polaritons, i.e. the well known normal-mode 
splitting [iSl [19] . As g increases, these peaks develop a 
marked asymmetry. For g ~ k, a clear splitting of the -I- 
peak occurs, corresponding to the resonant hybridization 
of one and two polariton states. Fig. [T] also shows results 
of a numerical (but non-perturbative) master-equation 
simulation [TJ] , showing our analytic approach is reliable 
even for moderately strong g. 

Large-detuned drives- The resonant-polariton inter- 
action is also interesting for drives far from the red- 
sideband, where the value of Gres "C wm- For a laser 
detuning near the minimum possible value A = —2ujm at 
which resonance is possible (and setting G = Gics), the 
polaritons are each either almost entirely phonon or pho- 
ton, implying a very small value of g ex gG/coM- However, 
as the — polariton is now almost purely phononic, its 
small damping rate and potentially large thermal occu- 
pancy enhances the self-energy in Eq. (11a) (i.e. k_ ~ 7, 
and n_ corresponds to the mechanical thermal occu- 
pancy) . We can quantify these effects by considering the 
value of pd[uj = E^], which will be suppressed by the 
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FIG. 2: Reflection coeflicient for a weak probe beam incident 
at a frequency u^-^ (defined in the lab frame), as measured 
in an OMIT experiment (upper inset). We take a one sided 
cavity with k^^/k, = 0.5, g = 0.5k. For each curve, A is 
labelled, and G — Grcs[A]. Remaining parameters are the 
same as Fig. 1. Lower inset: behaviour of g, Grcs and Ccfi as 
a function of detuning A of the main laser drive. 



hybdriziation physics described here. One finds: 



Pd[E+] = 



2/7T 



1 



K+ 1 + Ccff 



Ccff — 



4g2(l + 2n_) 



(14) 



For a large detuning, the effective cooperativity scales as 
Ceff ex C((7/wm)^, where C — AG'^/n'-f is the standard 
many-photon coupling cooperativity. Thus, in the large- 
detuned regime, resonant polaritons interactions allow 
one to amplify the effects of the nonlinear interaction by 
a factor C. Fig. [2] shows the evolution of the OMIT re- 
flection coefficient (which reflects the structure in Pd[w]) 
as the detuning A is varied, while keeping G tuned to 
the resonant value GreslA). 

Induced Kerr interaction- The nonlinear interaction 
in the resonant regime defined by Eqs. (I6])-(|7]) leads to 
a strongly enhanced two-particle interaction between — 
polaritons, mediated by the exchange of a + polari- 
ton (Fig. 3 ) . In a weakly-driven optomechanical sys- 
tem, Eq. ( 3 1 implies that phonons can mediate an effec- 
tive photon-photon interaction; however, as the virtual 
phonon is off-resonance, this interaction U ex g^ /ujyi- In 
contrast, the resonance condition E^ = 2E_ yields an 
induced interaction C/ros oc (j'^/k, an enhancement by a 
large factor c>c cjm/k. 

To assess the effects of the polariton-polariton in- 
teractions, we weakly drive our system with a second 
probe tone, and consider the 52 correlation functions 
g2u = {u^ u^ uu) / {u' u)"^ , where u = b,d,c+,C-. g2u is 
a measure of interaction induced correlations; (?2 < 1 sig- 
nifies non-classical correlation. Given the strong interac- 
tion experienced by — polaritons when the resonance con- 
dition E^ — 2E_is achieved, we expect that if the cavity 
is driven at the E— resonance, 32- will drop below 1. 
This is indeed the result of a numerical, master-equation 
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FIG. 3: Inset: Resonant interaction between — polaritons. 
Main: Numerically-calculated (?2 correlation function for — 
polaritons (<72-) and phonons {g^b), in the presence of an 
additional weak probe laser (frequency Ljpr)- Here, g — k, 
A = — ojm, G — 0.3tJM ~ Grcs. We have taken wm/k -^ 00 
to suppress non-resonant interaction effects. The probe am- 
plitude is e = 0.2k (52-), e ~ 0.3k {g2b)- Both phonon and 
polariton g2 functions drop below 1 due to the interactions, in- 
dicating non-classical correlations despite the fact g/uJM — 0. 
The dashed curve is the result of an analytic theory (see [14)V 



based calculation (see Fig. [s] and [Ilj). An analytic cal- 
culation based on a reduced state-space (similar to that 
in Ref. [T^) reproduces these results. For a weak probe 
drive at the E— frequency, it yields [13]: 



52- 



l + Ag'^/n^J 



(15) 



One also finds non-classical correlations for photons and 
phonons. Shown in Fig. |3]is the phonon 172 function g2b 
(for same parameters); it clearly drops below 1. The 
double-peak structure of this curve is the result of the 
drive inducing correlations between — and -|- polaritons; 
it also occurs in the behaviour of {wb) (see EPAPS for 
more details [T4]). 

Conclusions- We have presented a systematic ap- 
proach for describing nonlinear interaction effects in 
a driven optomechanical system, identifying a regime 
where a resonance enhances interactions between polari- 
tons. We have discussed how this would manifest itself 
in a OMIT-style experiment, as well as in (72 correlation 
functions. The polariton interactions we describe could 
be extremely interesting when now considered in lattice 
systems, or when considering the propagation of pulses. 

We thank W. Chen and A. Nunncnkamp for useful dis- 
cussions. This work was supported by CIFAR, NSERC 
and the DARPA ORCHID program under a grant from 
AFOSR. Note added- During the preparation of this pa- 
per, we became aware of a related work by B0rkje, Nun- 
ncnkamp, Teufel and Girvin. 
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SUPPLEMENTAL INFORMATION 

LANGEVIN EQUATIONS FOR LINEARIZED OPTOMECHANICAL SYSTEM IN THE POLARITON 

BASIS 



We start by diagonalizing the linearized optomechanical Haniiltonian given in Eq. (2) of the main text, working as 
always in a displaced interaction picture set by the laser drive on the cavity. Introducing 



X :^ [b d b^ S ] , Y ^[c^ c+ c 
the diagonalization can be expressed in terms of a 4 x 4 real matrix U: 

Y = V -X. 






(SI) 



(S2) 



U can be found by standard means, though its general form is both cumbersome and unenlightening; we define 
V = U~^. It is slightly less unwieldy in the special (but relevant) case of a drive at the red-detuned mechanical 
sideband A = — wm- In this case, we have simply: 



1 



c± 



while the inverse transformation is defined by: 

V2j 



{E± - LOM){d^ ±b^) + {E±+ LOM){d±b) 



a/8cJm 
1 



a/8cJm 






(S3) 

(S4) 
(S5) 



Next, we include the coupling of our system to the cavity and mechanical dissipative baths in the standard way, 
treating these baths as Markovian over frequencies of interest. Consider first the coupling to the cavity bath. Prior 
to making displacement and interaction-picture transformations, the system-bath coupling has the form (see, e.g., 
Ref.H]): 



H^^Y.^J]I,, H,,nt = -*y^E(/^-/i)(^ + «^) 



(S6) 



where fj is the lowering operator for a bath mode, and pc is the density of states of bath modes (which we treat 
to be constant over frequencies of interest). In now moving to an interaction picture at the laser frequency, we 
also transform the bath modes, i.e. fj — > /je^"^^*. Formally, the interaction picture transformation involves a 



unitary U{t) = exp 



-iuJLt 



;t; 



\a' a + J2i fj fj) ^nd transforms the bath Hamiltonian to iJ„ = J2ji^j ^ '^L)f ] fj- 

In this interaction picture, the counter-rotating terms in H^i„t will explicitly oscillate at ±2^^ = ±2(wc + A). 
Even if we now write our photon operators a in terms of polariton operators, there is no possibility of having these 
terms becoming resonant, as the cavity frequency is much larger than any other frequency scale in the problem (i.e. 
Wc ^ \A\,ijJm, E^, E-). As such, one can safely make a rotating- wave approximation in the photon basis, resulting 
in a standard system-bath interaction which is stationary in the interaction picture: 



H.- 



^EC/I"'--^*/. 



27rp, 



^ Y. [j'l (^2lC- + V22C+ + V23CI + V24cl) - h.C 



2Trpc 



(S7) 



We have also made the displacement tran sformation a — d + a as discussed in the main text. Note that in writing d in 



terms of polariton operators via Eq. (S2), we obtain terms of the form /„-cJ., which can cause polariton heating even 



t;>t 



if the cavity bath is at zero temperature. Such terms are physical and must be retained. Formally, In the interaction 
picture the bath now has negative frequency modes which can make such processes resonant. In more physical terms. 



the combination of zero-point bath fluctuations with the cavity driving can excite the polaritons. This mechanism 
has been discussed in other contexts under the name "quantum activation" by various authors [21 [3] . 

We turn now to the mechanical bath, where the basic interaction Hamiltonian can be written in an analogous way: 



H-,--J2'^39]9v 



H. 



-y.int 



7 



27rp„ 



E(^.-^l)(^ 



&t 



(S8) 



Here, ijj is a lowering operator for a mechanical bath oscillator, and Pm is the density of states of mechanical bath 
modes (which we also treat to be constant over frequencies of interest). As there is no direct driving of the mechanical 
resonator, the analysis here is somewhat simpler. We first re-write the phonon operator b in the polariton basis, and 
then make a rotating-wave approximation. The justification is that counter-rotating terms such as g]cl can never be 
made resonant; as the mechanical resonator is not driven, there is no quantum activation mechanism involving the 
mechanical bath. We thus obtain 



H, 



7,m£ 



e/5^i:(s]-+H 



where (V = U^i) 



«:5=7(14i + ^3i) , 4' = 7 (^12 + V^32) 



(S9) 



(SIO) 



Having now established the correct system-bath coupling Hamiltonians in the polariton basis we wish to use, we 
can derive the Hcisenberg-Langevin equations for our system in the standard manner. As each bath couples to both 
+ and — polaritons, one obtains off-diagonal damping terms, (e.g. the mechanical bath produces a force on the -|- 
polariton that is proportional to the amplitude of the — polariton). Such terms will dynamically couple + and — 
polaritons, and can be included in our theory in a straightforward manner (i.e. by including Green functions that are 
off-diagonal in the polariton index). However, in the weak-dissipation limit of interest {k,j ^ E^,E_, \E^ — E_\), 



the mixing effects induced by such terms is strongly suppressed, 
resulting in the form: 



As such, we will drop off-diagonal damping terms. 



dtCait) = 
where the polariton damping rates Ko- are given by: 



K^'+K 



(T^2l)' - {V23f 



iEa + -^j Co-(t) - 0vCcr(i), 



kI' +K 



(1^22)' - (1^24)' 



(Sll) 



(S12) 



The noise operators ^(j (t) are each linear combinations of the input noise emanating from the mechanical and cavity 
baths. In the interaction picture we use, we will be sensitive to noise in the cavity bath at frequencies near Wc, 
and noise in the mechanical bath at frequencies near ujm- In the limit of interest where the physical temperature 
T ^ hwc/kB, there will be no thermal noise incident from the cavity bath at the frequencies of interest. Also, as 
we focus on regimes where the polariton damping rates are much smaller than their energies, we can treat the noise 
operators as being white noise (as is standard in input-output theory treatments). We thus have 



{il{t)Ut')) = nj [t - t') , {Ut)il{t')) = (1 + ".) 5 (t - t') , 
where the effective temperatures of the two noises are given by: 

n_ = -^(K*^nB[S_] + (y23)'K), fl+ = -^ (k^ ns [^+] + (V24)' k) 



(S13) 
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ns [Ea\ denotes the Bose-Einstein distribution function evaluated at energy E„ and the mechanical bath temperature. 
Finally, one also finds that ^+ and ^_ are correlated with one another. Similar to the situation of off-diagonal 
damping terms, such noise correlations could easily be included in our theory; however, as they play no role in the 
regime of interest where -Eo-, \E+ — E^\ 3> ^,7, we drop them in what follows. One also finds that anomalous noise 
correlators can be non-zero (e.g. (C(t(^)?o-'(0)))- Again, while such terms can be retained in our theory, they play no 
role in the weak-damping regime of interest, and we hence drop them in what follows. 



UNPERTURBED POLARITONS GFS 



The standard definitions of the three relevant GFs needed in the Keldysh technique are (ct stands for ±, but similar 
definitions hold for c^ = 6 or d): 



G^„, M = -I I dte{t)l c,{t),cl,{Q) 

-oo 

dte{-t)i\c„{t),cl,{Q) 

Gf,,H ^ -I /"dt({c,(i),4,(0)})e 



Gt' H 
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(S16) 
(S17) 



The retarded and advanced GFs keep track of spectral information, whereas the Keldysh Green function G^ also 
keeps track of the occupancy of states. For the linearized (non-interacting) theory, the GFs are easily obtained from 
the Langevin equations in Eq. (Sll ). These free GFs (which we denote by Q) are diagonal (only non-zero for a = a') 
and given by: 



G^A^] = 



uj- E„ +i- 



g'^A^] - (2n. + i)(t;f,H-e..H), 



(S18) 
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Finally, as there are no off-diagonal polariton Green functions, we can use Eq. (S2) to write the photon retarded 
Green function G^^[a;] as: 



GUA = E {G'tG^^M + D, [G^,[- 



(S20) 



where 



G^ 



V^\, C+^Vl, D^^Vi,, D+ = V, 



(S21) 



As the Green functions remain diagonal in the polariton index even with the nonlinear interaction (in the regime of 
interest, see below), the above relation also holds for the full Green functions (i.e. including the self-energy associated 
with g). 



H2 IN THE POLARITON BASIS 



Using the change of basis matrix V = U ^ (c.f. Eq.(S2)), we can re- write the non-linear interaction H2 in Eq. (3) 
of the main text in the polariton basis via: 



H2=9 
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A,B 



Expanding this equation allows one to obtain the interaction coefficients g^^,g.,, in terms of g and matrix elements of 
V. In particular, the coefficient g of the resonant interaction process c_c_c+ will be given by: 



~9 = 5(1^22 V^21 (^11 + ^13) + ^^23^24(^11 + Vu) + y23"^'2lC^14 + ^^12)). 



(S23) 



As the normal-mode transformation described by U mixes raising and lowering operators, H2 will not be normal- 
ordered when written in terms of polariton operators (even though it is normal ordered when written in terms of 
photon and phonon operators). Normal-ordering H2 in the polariton basis yields the form: 



H2 =: H2:+ (A_c_ + A+c+ + h.c.) 



(S24) 



where the colons indicate normal-ordering in the polariton basis, and the constants A^ oc g. The first term is the 
normal-ordered polariton interaction written in Eq. (5) in the main text. We next make a unitary displacement 



OJ + UJi 



a) E^_fa;l = i 



E« M 




b) S«+H = 




> ,t> >:t> 

77X^1 



fTq dpi! dil 



e-M 



^++M 



FIG. 4: (a) Feynman diagrams describing the dominant + and — polariton self-energies (to second order in g) in the case 
when the resonance condition E+ — 2_B_ is met. We have shown diagrams for the retarded self-energies; similar diagrams 
also determine the advanced and Keldysh self-energies. The structure of propagators in Keldysh space is indicated by writing 
explicit classical (d) and quantum (g) indices at the ends of each Green function, (b) The non-resonant interaction terms induce 
effective two-particle interactions between polaritons, as described by -ffjvu in Eq. (8) of the main text. These two-particle 
interactions modify the resonant self-energies depicted in Fig. Q; shown above are the diagrams describing the modification 

of the -I- polariton self-energy by th e effec tive interaction U The summation of these ladder diagrams results in the simple 

energy shift indicated in Eqs. (S27l,(S28l. 



transformation of the form Ca -^ Cg- + £„ 
Hamiltonian; to leading order in g, c^ ^ 



, where the constants c^ are chosen to eUminate all linear-in-CCT terms in the 
—Aa-/Ea-- The resulting coherent Hamiltonian has the form: 



Ho 



Hn 



st At 



h.c. 



Here Zi, Z2 are 4x4 matrices whose entries are all order g'^/Eg- ^ g'^/wM- The quadratic terms on the RHS describe 
corrects to the linear Hamiltonian Hi — ^^ Ecc\ca arising from the displacement transformation. In principle, one 
could combine these with the terms in iJi, re-diagonalize the resulting quadratic Hamiltonian, obtaining a new basis 
of non-interacting polaritons. However, to leading order in g^/wM, all that is needed is to retain the diagonal elements 
of Zi, which simply shift the polariton energies but do not change their wavef unctions. These energy shifts (denotes 
e+,e_) are just absorbed into the definition of our free Green functions. We are then left with a normal-ordered 
interaction Hamiltonian H2 in the polariton basis which can be addressed perturbatively. 



PERTURBATIVE TREATMENT 



We have calculated the full Keldysh self-energy corresponding to H2 to order g'^ , without any further approximation. 
We have also done this calculation in the case where G is not so large, such that one should work in the original basis 
of photons and phonons. Results of these full calculations will be presented elsewhere. Here, like in the main text, 
we will focus on the resonant-interaction regime described in the main text, where the condition Ej^ — 2E_ enhances 
certain scattering processes. In this regime, the dominant self-energy processes for -I- and — polaritons are depicted 
in Fig. H Using the standard rules of the Keldysh technique [3] , these diagrams correspond to: 



2,52 / ^ ^gK_ [^,] gn 



27r 



u} - (E+ - E^) + I — ^-^ 



s^+ M 



2^g^ I '^g-_ W] g- [. - .'] = 2g^ ^^g;^ 

' ZTT U) — 2E_ + IK_ 



(S25) 
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Note that because of the resonance condition, these self-energies scale as g'^/n, whereas all other self-energy diagrams 
are suppressed by an additional small parameter k/ojm- They can thus be neglected in the limit k/ujm — >■ 0. For a 
small but realistic value of k/ujm, the small-energy shifts associated with the non-resonant interaction terms in H2 
can shift the condition needed for resonance. To describe these small shifts (which can be important for realistic 
parameters), it is sufficient to use standard second-order perturbation theory to treat the non-resonant terms. This 
is conveniently done via a Schrieffer- Wolff transformation, where a unitary transformation is performed to eliminate 
the non-resonant terms to leading order in the Hamiltonian. The procedure generates new terms however at second 
order in g. Keeping only such terms which do not change the total number of polaritons, we obtain the general form 
given in the Hamiltonian Hjyjf (Eq. (8) of the main text). Note that in this equation, the energy shifts S^r include 
both the shifts arising from the non-resonant terms, as well as the energy shifts e^ coming from the normal-ordering 
procedure. Both such terms scale as g'^/ujM- 
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Including the energy shifts associated with the non-resonant processes (as described by Hnb. in Eq. (8) of the main 
text), the self-energies above are modified to: 



4g^ 



K_-|-K I 



2g 



Lo — E^ + i 

2 (l + 2n_) 
ijj — E + iK_ 



E+- =E+-E^+5+-8^+ U+- (n_ - n+) , 



i^__ =2{E^+ 5_) + 2U— (1 + 2n_] 



(S27) 

(S28) 



The corrections due to the Kerr- type interaction constants U^ and U can be obtained by including ladder diagrams 

in the self-energies, as shown in Fig. [4J3. 

OMIT REFLECTION COEFFICIENT 

In an OMIT style-experiment, in addition to the main driving laser (which gives rise to the many-photon interaction 
G), a second weak drive tone (the "probe") is applied at a frequency Wpj- to the cavity. This driving is described by 
a term in the Hamiltonian: 



i/n 



' ^cp I ^ ^'u 



3~'"P'* -f h.C. 



(S29) 



where we work in the interaction picture determined by the main drive laser frequency, and hence Wpi- — Wpj. — Wi . 
Kcp parameterizes the coupling of the drive port to the cavity: for a one-port cavity, the total cavity damping rate 
'i — i^c-p + '*int, where Kjnt describes internal cavity losses. 

As the amplitude din probe drive is weak, we can use standard linear response theory to describe its effects. The 
Kubo formula thus tells us that to first order in din.pr, the change in the cavity amplitude will be given by: 



5 id{t)) = -i^Ti^ r dt' (Ji„,p.G« (t - i')e-^"-*' - dl^^.G%{t - t')e^"-*') 



— OO 

— iu)-r,rt 



-zV^(e-*"-*di„,p,Gf,[(ip,] ~ e--*dl,p,G«-hc:;p,]) 
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Here, G£j[a;] is the retarded Green function of the cavity as defined above, calculated to zeroth order in ffpr, but 
including the effects of the nonlinear interaction g. G^^\io\ is the corresponding anomalous Green function defined as: 



(S32) 



GIj[lo] = -ij dt6{t){^\d{t),d{Q) 
Now, the standard input-output relation between input, output and cavity fields is [T]: 

dont{t) = dinit) + y/K^d{t) 



(S33) 



Taking the average value of this equation, and defining the elastic amplitude reflection coefficient r[a;pr] as the 
amplitude of ((iout(i)) at the probe frequency divided by din, we obtain 



as given in the main text. 



I - iKcpG^di^pr], 



MASTER EQUATION SIMULATION 



(S34) 



Starting from the system-bath Hamiltonians written in the polariton basis (Eqs. (S7) and (S9|), one can trace over 



the dissipative baths and derive a master equation for the reduced density matrix p describing the polaritons using 
standard Born-Markov approximations [S]. One obtains: 



dtp=-i Hi + H2,p + Cp 



(S35) 
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where the coherent system Hamiltonian Hi + H2 is written without any approximation, and the super-operator C 
describes the effects of the dissipative baths via standard Lindblad terms, 



C = K_(l + n_)L[c_] + K_n_L[cl] + k+(1 + n+)L[c+] + K+n+i[c^] 



(S36) 



with 



m- 



= c ■ c* 



-Ac^c,-}. 



(S37) 



Note that this master equation corresponds to each polariton seeing independent thermal baths, in direct analogy to 



the form of the quantum Langevin equations in Eq. (Sill 



While Eq. ( S35 ) is not a convenient starting point for deriving analytic results, it does allow us to numerically study 



the system without having to assume a small value of g. Using Eq. (S35) and the quantum regression theorem [S], 
we have numerically calculated the cavity density of states Pd[w] = — Im G^^[w]/7r, finding good agreement with our 
analytic perturbative results even for g as large as k (see Fig. 1 in main text). Note that to find agreement with these 
numerical results, it was crucial to i nclud e in t he a nalvtic theory the corrections associated with the non-resonant 
interaction processes, Hfqp; (c.f. Eqs.(S27) and (S28)). 



Finally, we have also numerically studied our system using a more conventional master equation, in which the 
dissipation is described by Lindblad superoperators which act in the photon and phonon basis, i.e. replace C in 
Eq. (S35) with Cq, where: 



£0 = KL\d] + 7(1 + nM)L[h] + -ffiMLp], 



(S38) 



and fiM is a Bose-Einstein distribution evaluated at the mechanical frequency wm and mechanical bath temperature. 
For the parameters studied in the paper, this conventional master equation yields results very similar to those obtained 
from Eq. 



(S35I 



32 CORRELATION FUNCTIONS 



As discussed in the main text, the resonant two-particle interaction between — polaritons can lead to non-classical 
values of the 52 correlation functions. In Fig. 3 of the main text, we have used a numerical solution of our master 
equation and the quantum regression theorem to calculate 32 functions of both photons and — polaritons, taking 
f; = K, A = — ojm and tuning G to Gics — 0.3wm to ensure the resonance condition E^ — 2£'_. For simplicity, we also 
take the large wm limit (i.e. lum/i'^ -^ 00, g/^M ^ 0), so that corrections due to non-resonant interaction terms in H2 
can be neglected. We have also included a weak probe drive on the cavity at a frequency cjpr close to E- by including 



a term of the form H^^ (c.f. Eq.(S29 1). When writing this driving term in terms of polarities, we will have terms that 



effectively drive both — and + polaritons. However, in the large wm limit and for Wpi- ~ E^, the direct driving term 
on the -|- polaritons will be strongly off-resonance and can be neglected. The probe field driving Hamiltonian thus 
reduces to: 



i/pr = ee 



A^l 



(S39) 



Eq. (15) of the main text and the dashed curve in Fig. 3 are the results of a simple analytic theory which accurately 
describes the — polariton (72 function in the limit of an extremely weak probe drive amplitude e. For a weak drive 
(and only keeping the resonant interaction process in H2), we can restrict attention to the four lowest eigenstates of 
the coherent Hamiltonian. These are: 



|V'o) = |0,0), 

IV'i) = |i,o), 



^0 = 0, 
El = £:_ 

E^^U 



\iP2) = - sin 6l|0,l)+ cos 6*12,0), 
|i^3)-cos0|O,l)+sin0|2,O), E3 = h^E^ 



E^ 
E^ 



{\E, 



E^] 



2.92, 



C^E+-E^)^ + 2~g^, 
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polarito ns, a nd tan 20 = 2\/2g/{E+ — 2E- 



where |n_,n_|_) denotes a Fock state having 7i_ — polaritons and n^ 

The system is driven by adding the probe drive Hamiltonian Hp,- given in Eq. (S39). Including the drive, the wave 

function is written in terms of the above states as: 



1^) = Aoo |0, 0) + Aio |1, 0) + A20 |2, 0) -I- ^01 10, 1). 



(S44) 
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FIG. 5: Plot of average polariton, photon and phonon numbers versus probe frequency ujp-c, for para meter s identical to those 
used in Fig. 3 of the main text. We have taken a drive strength of e = 0.2fi;, where e is defined in Eq. (S39l. The double-peak 



structure in the phonon population is due to higher-order effects of the probe drive which induce correlations between -I- and 
— polarities. 



We take into account the dissipation with the anti-Hermitian Hamiltonian H^g^ 



— — ii(K_c_c_ + nj^c}j_Cj^ 



the rotating frame corresponding to the unitary operator e^'-'^-'^^'^'^'^+'^+'^p'^*, the Hamiltonian reads 
iJRWA = (A- - \n-)c^_c- + (A+ - \k+)c^_^c+ + g{clc^c^ + c+clcl) + e{cl + c_), 



In 



(S45) 



with the detuning A_ = E_ — LUp,. and A+ = E^ — 2ujpi. = 2A_ at the resonance E^ — 2_E_. The dynamics is 
given by dt\ip) = —iH\'ip). Ignoring other states than those involved in j?/;), the steady state number and correlation 
(72- = (cLclc_c_)/(c_c_)^ are 



{clc\ 



92- 



A1 + (k_/2)2' 

[A^ + (At_/2)2][4A^ + (Ac_/2)2] 

[2A2_ - 52 _ («;_/2)2]2 + 9A2_(k_/2)2 ' 
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where we approximate Aqo = 1, take k+ = k_ and keep the leading terms in e. The minimum value of (72- is found 



at Wpr = E^ 



with (72- — n't/i'ig'^ + ^2^)2 and the correlation tends toward unity for infinite detuning. 



We comment briefly on the unusual double-peak resonance exhibited by 52b ['^pr] in Fig. 3 of the main text. This 
structure is also present in the behaviour of the average phonon number versus probe drive but not in the behaviour of 
the average photon number, see Fig. [5] This difference is a higher-order effect of the probe drive-field that would not 
be captured in linear-response (i.e. contributions to (b^b) beyond order e^). At higher orders, the probe field causes 



non-zero correlations between -I- and — polaritons, e.g. {c\_C-) ^ 0. It follows from the normal-mode transformation 



of Eq. (S2| that such averages contribute to the photon and phonon populations with opposite sign; this ultimately 



explains the strong difference in the behaviour of photon and phonon number versus u>- 



pr- 
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